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Abstract 

The paper investigates the properties of certain biorthogonal polynomials appearing in a specific 
simultaneous Hermite-Pade approximation scheme. Associated to any totally positive kernel and a 
pair of positive measures on the positive axis we define biorthogonal polynomials and prove that 
their zeroes are simple and positive. We then specialize the kernel to the Cauchy kernel and 
show that the ensuing biorthogonal polynomials solve a four-term recurrence relation, have relevant 
Christoffel-Darboux generalized formula and their zeroes are interlaced. In addition, these poly- 
nomial solve a combination of Hermite-Pade approximation problems to a Nikishin system of order 
2. The motivation arises from two distant areas; on one side, in the study of the inverse spectral 
problem for the peakon solution of the Degasperis-Procesi equation; on the other side, from a random 
matrix model involving two positive definite random Hermitian matrices. Finally, we show how to 
characterize these polynomials in term of a Riemann-Hilbert problem. 
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7 Acknowledgments 

A Appendix: Proof of Extended Christoffel-Darboux Identities 

1 Introduction and motivations 

This paper mainly deals with a class of biorthogonal polynomials {pn{x)}n, {9n(y)}N of degree n satisfying 
the biorthogonality relations 

/ / Pn[x)qm{y) ■ =Omn, (1-1) 

where da, d/3 are positive measures supported on M_|- with finite bimoments. These polynomials will be 
introduced in Sec. [2]in a more general context of polynomials associated to general totally positive kernels 
(Def. 12. 1[) with which they share some general properties in regard to their zeroes. 

While these properties are interesting in their own right, we wish to put the work in a more general 
context and explain the two main motivations behind it. They fall within two different and rather distant 
areas of mathematics : peakon solutions to nonlinear PDEs and Random Matrix theory. 



Peakons for the Degasperis-Procesi equation. In the early 1990's, Camassa and Holm [TT] in- 
troduced the (CH) equation to model (weakly) dispersive shallow wave propagation. More generally, the 
CH equation belongs to the so-called b-family of PDEs 

ut - Uxxt + ib+ l)uuj. = buxUxx + uuxxx, (x, i) e M^, & e R, (1-2) 

Two cases, 6 = 2 and 6 = 3 within this family are now known to be integrable: the case 6 = 2 is the 
original CH equation whereas the case 6 = 3 is the Degasperis-Procesi [14] (DP) equation, which is more 
directly related to the present paper. 
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In all cases the b-family admits weak (distributional) solutions of the form: 

n 

«(x,t)=^m,(Oe-l----(*)l, (1-3) 
if and only if the positions Xi{t) and the heights mi{t) satisfy the system of nonhnear ODEs: 

n n 

ifc = ^m,e"l'^'="'''l, riifc = (6- l)^TOfcTOiSgn(a;fc -a::,)e~l^''--^*l, (1-4) 

i=l i=l 

for fc = 1, . . . , 71. The non-smooth character of the solution manifests itself by the presence of sharp peaks 
at {xk\, hence the name peakons. For the CH equation the peakons solution were studied in [51[T], while 
for the DP equation in PHlllIj; in both cases the solution is related to the isospectral evolution of an 
associated linear boundary-value problem 



b = 2 [CH) 


5 = 3 {DP) 




-0"'(e,z) = z.9(O0(e,^) 


(^(-1) = 0(1) = 


0(-l) = 0'(-l) = (/.(I) = 



(1-5) 



The variables x, ^ and the quantities m, g, u are related by 

C = tanh^^^^ , g(C) ^ ^ ^ m{x) , m{x,t) = u{x,t) - u^^{x,t) (1-6) 

Because of the similarity to the equation of an inhomogencous classical string (after a separation of 
variables) we refer to the two linear ODEs as the quadratic and cubic string, respectively. The case of 
peakons corresponds to the choice 

n 

m{x,t) = 2^S{x - x^{t))m,(t) (1-7) 

The remarkable fact is that in both cases the associated spectral problems have a finite positive spectrum; 
this is not so surprising in the case of the quadratic string which is a self-adjoint problem, but it is quite 
unexpected for the cubic string, since the problem is not self-adjoint and there is no a priori reason for 
the spectrum to even be real [H] . 

As it is natural within the Lax approach to integrable PDEs, the spectral map linearizes the evolution 
of the isospectral evolution: if {zj} are the eigenvalues of the respective boundary value problems and 
one introduces the appropriate spectral residues 

^.-res^d., Wiz):=^ (1-8) 



Z — Zj 
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then one can show [20j that the evolution hnearizes as foUows (with the dot representing the time 
evolution) 

ife - , ^ = 1 (1-9) 

Since this is not the main focus of the paper, we are deliberately glossing over several interesting points; 
the interested reader is referred to [21] and our recent work [S] for further details. In short, the solution 
method for the DP equation can by illustrated by the diagram 



I DP flow I 



evolution of the extended spectral data 
r /,\ /j-\\n inverse spectral map {Zfc(t)=Zfc 

In the inverse spectral map resides the role of the biorthogonal polynomials to be studied here, as we 
briefly sketch below. The inverse problem for the ordinary string with finitely many point masses is 
solved by the method of continued fractions of Stieltjes' type as was pointed out by M.G. Krein ([17]). 
The inverse problem for the cubic string with finitely many masses is solved with the help of the following 
simultaneous Hermite-Pade type approximation (\21\) 

Definition 1.1 (Pade-like approximation problem). Let dfi(x) denote the spectral measure associated 
with the cubic string boundary value problem and l^Isi = J -L-.d^[x), §■ = J ^;^d^{x)d^{y) denote 
the Weyl functions introduced in \21^ . Then, given an integer I < k < n, we seek three polynomials 
{Q, P, P) of degree k — \ satisfying the following conditions: 

[Approximation] :W=^ + (^-^) , ^ = ^ + ^ (z ^ oo). 

[Symmetry].- Z* Q + W* P ^ P = O {z ^ oo) withW*{z) = -W{-z), Z*{z) = Z{~z). 

[Normalization]; P(0) = 1, P(0) = 0. 

This approximation problem has a unique solution ([21]) which, in turn, is used to solve the inverse 
problem for the cubic string. We point out that it is here in this approximation problem that the Cauchy 
kernel ^^^^ makes its, somewhat unexpected, appearance through the spectral representation of the second 
Weyl function. 

Random Matrix Theory The other source of our interest in biorthogonal polynomials comes from 
random matrix theory. It is well known [22| that the Hermitean matrix model is intimately related to 
(in fact, solved by) orthogonal polynomials (OPs). Not so much is known about the role of biorthogonal 
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polynomials (BOPs). However, certain biorthogonal polynomials somewhat similar to the ones in the 
present paper appear prominently in the analysis of "the" two-matrix model after reduction to the 
spectrum of eigenvalues O El O I15| ; in that case the pairing is of the form 

Pnix)q,niy)e'''^da{x)dP{y) = Smn, (1-10) 

and the associated biorthogonal polynomials are sometimes called the Itzykson-Zuber BOPs, in short, 
the IZBOPs. 

Several algebraic structural properties of these polynomials and their recurrence relation (both mul- 
tiplicative and differential) have been thoroughly analyzed in the previously cited papers for densities of 
the form da{x) ~ e^^'-^^^dx, d/3(j/) = e^^^^^My for polynomials potentials Vi{x), V2{y) and for potentials 
with rational derivative (and hard-edges) in [3]. 

We recall that while ordinary OPs satisfy a multiplicative three-term recurrence relation, the BOPs 
defined by (|l-10p solve a longer recurrence relation of length related to the degree of the differential 
dVj{x) over the Riemann sphere [3]; a direct (although not immediate) consequence of the finiteness 
of the recurrence relation is the fact that these BOPs (and certain integral transforms of them) are 
characterized by a Ricmann-Hilbert problem for a matrix of size equal to the length of the recurrence 
relation (minus one). The BOPs introduced in this paper share all these features, although in some 
respects they are closer to the ordinary orthogonal polynomials than to the IZBOPs. 

The relevant two-matrix model our polynomials arc related to was introduced in |10| . We now give a 

(2) 

brief summary of that work. Consider the set of pairs Tl\_ := {(Mi, M2)} of Hermitean positive- definite 
matrices endowed with the (C/(iV)-invariant) Lebesgue measure denoted by dMidM2. Define then the 
probability measure on this space by the formula: 

dM(Mi,M2) = ^ det(Mi+M2)^ ('"^'^ 

(2) 

where Zj^ (the partition function) is a normalization constant, while a'(Mi), /3'(M2) stand for the 
product of the densities a' , P' (the Radon-Nikodym derivatives of the measures da,d/3 with respect to 
the Lebesgue measure) over the (positive) eigenvalues of Mj. 

This probability space is similar to the two-matrix model discussed briefly above for which the coupling 
between matrices is e^"^'^^!*^" [H] instead of det(7\/i + Ah)^^ ■ The connection with our BOPs dTT]) 
is analogous to the connection between ordinary orthogonal polynomials and the Hermitean Random 
matrix model whose probability space is the set of Hermitean matrices Hn equipped with the 
measure d/ii(M) := —jjja'{M)dM. In particular, we show in [10] how the statistics of the eigenvalues 
of the two matrices Mj can be described in terms of the biorthogonal polynomials we are introducing 
in the present work. A prominent role in the description of that statistics is played by the generalized 
Christoffel-Darboux identities we develop in Section |31 

We now summarize the main results of the paper: 
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- for an arbitrary totally positive kernel K{x,y) and arbitrary positive measures da,d/? on we 
prove that the matrix of bimoments lab ■= jrig2 x°'y^K{x, y)da{x)df3{y) is totally positive (Thm. 

HOD; 

- this implies that there exist, unique, sequences of monic polynomials of degree n, p„(a;),g„(y) 
biorthogonal to each other as in (|2.ip : we prove that they have positive and simple zeroes 
(Thm. [2751) : 

- we then specialize to the kernel K{x,y) = this case the zeroes of Pn{x) (9n(y)) are 
interlaced with the zeroes of the neighboring polynomials (Thm. 13.21 ); 

- they solve a four— term recurrence relation as specified after (|l-ip (Cor. [472]); 

- they satisfy Christoffel Darboux identities (Prop. HS], Cor. \4M Thms. [HH [53)) 

- they solve a Hermite-Pade approximation problem to a novel type of Nikishin systems (Sec. [5l 
Thms. [5lll[5l2|); 

- they can be characterized by a 3 x 3 Riemann— Hilbert problems, (Props. 16.11 16. 2p ; 

In the follow-up paper we will explain the relation of the asymptotics of the BOPs introduced in this 
paper with a rigorous asymptotic analysis for continuous (varying) measures da, d/3 using the nonlinear 
steepest descent method [9]. 

2 Biorthogonal polynomials associated to a totally positive ker- 
nel 

As one can see from the last section the kernel K{x,y) = -^^,x,y > 0, which we will refer to as the 
Cauchy kernel, plays a significant, albeit mysterious, role. We now turn to explaining the role of this 
kernel. We recall, following [12], the definition of the totally positive kernel. 

Definition 2.1. A real function K{x,y) of two variables ranging over linearly ordered sets X and y, 
respectively, is said to be totally positive (TP) if for all 

Xi < X2 < ■ ■ ■ < Xm, yi < y2 < ■ ■ ■ < Urn € A", ?/j G 3^, m £ N (2-1) 

we have 

det[if(a;„%)]j<,,^.<,„ >0 (2-2) 
We will also use a discrete version of the same concept. 
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Definition 2.2. A matrix A :— [aij], i, j = 0, 1, • • - n is said to be totally positive (TP) if all its minors 
are strictly positive. A matrix A := [o-ij], i, j = 0, 1, • • • n is said to be totally nonnegative (TN) if all its 
m,inors are nonnegative. A TN matrix A is said to be oscillatory if some positive integer power of A is 
TP. 

Since we will be working with matrices of infinite size we introduce a concept of the principal trunca- 



Definition 2.3. A finite n + \ by n + 1 matrix B := \bi,j\, h j = 0, 1, • • • n is said to be the principal 
truncation of an infinite matrix A :~ [a^], i, j = 0, 1, • • • if bij = aij, i,j — 0,1, - ■ ■ n. In such a case B 
will be denoted A[n\. 

Finally, 

Definition 2.4. An infinite matrix A := [a^], i,j ~ 0, 1, • • • ts said to be TP (TN) tf A[n] is TP (TN) 
for every n = 0, 1, • • • . 

Definition 2.5. Basic Setup 

Let K{x,y) be a totally positive kernel on x and let da,d/3 be two Stieltjes measures on 
R+. We make two simplifying assumptions to avoid degenerate cases: 

1. is not an atom of either of the measures (i.e. {0} has zero measure). 

2. a and (3 have infinitely many points of increase. 
We furthermore assume: 

3. the polynomials are dense in the corresponding Hilbert spaces Ha := L^(M_|_, da), Hp :~ i^(R+, d/3). 



4. the map K : Hp Ha, Kq{x) :~ j K {x , y)q{y)dP{y) is bounded, injective and has a dense range 
in Ha . 

Under these assumptions K provides a non-degenerate pairing between Hp and Ha- 



Remark 2.1. Assumptions 3 and 4 could be weakened, especially the density assumption, but we believe 
the last two assumptions are the most natural to work with in the Hilbert space set-up of the theory. 

Now, let us consider the matrix 2 of generalized bimoments 



tion. 




(2-3) 




(2-4) 



Theorem 2.1. The semiinfinite matrix 2 is TP. 
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Proof. According to a theorem of Fckete, (see Chapter 2, Theorem 3.3 in [TH] ), we only need to consider 
minors of consecutive rows/columns. Writing out the determinant, 

^n' '^Ct[Ia+i,b+j]o<i,j<n-l 

we find 

= // \{^';y]\{^r'yr'K{x„y,)A''c^iX)d-p{Y)^ 

<tGS„ j=l j=l 

« » n n 

// C{xrC(Y)''A{X) n yf' n ^^(^.>%)d"ad"/3. 

Since our intervals are subsets of M+ we can absorb the powers of C{X),C{Y) into the measures to 
simplify the notation. Moreover, the function S{X,Y) := YVj=i ^i^j^yj) enjoys the following simple 
property 

S{X,Y,) = S{X,-.,Y) 

for any a G Sn- Finally, the product measures d"a = d"a(X),d"/3 = d"/3(y) are clearly permutation 
invariant. 

Thus, without any loss of generality, we only need to show that 



A. JJ A(X)[]yf i5(X,r)d"ad"/3 > 0, 



i=i 

which is tantamount to showing positivity for a = b ~ 0. First, we symmetrize D„ with respect to the 
variables X; this produces 



^» = ^ E JJ A(X.)n2;f '^(X„,y)d"ad"/3- 1 ^ JJ A{X)eia)l[yf'S{X,Y,-^)d-ad-p = 
-,Y. II A(X)e(a)f[yi-i5(X,y)d"ad"/3 = iy jj A(X)A(r)5(X, r)d"ad"/?. 



n! ^-^ II LL j TV. 1 1 

Subsequent symmetrization over the Y variables does not change the value of the integral and we obtain 
(after restoring the definition of S{X,Y)) 

'-"■^ Tes„ j=i 

JJ A(X)A(y)det[/^(x„yj)],,,<„d"ad"/?. 

Finally, since A(X)A(y) det[if (x^, j/j)]ij<„d"Q;d"/3 is permutation invariant, it suffices to integrate over 
the region < xi < X2 < • • • < a;„ x < j/i < ?/2 < • • • < and, as a result 

Dn= li A(X)A(r)det[/f(a;„yj)],,,<„d"ad"/3. (2-5) 

0<yi<y2<-<Vn 
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Due to the total positivity of the kernel K{x, y) the integrand is a positive function of all variables and 
so the integral must be strictly positive. □ 



To simplify future computations we define [a;] := (1, 



so that the matrix of generalized 



bimoments (|2-4p is simply given by: X = {[x 



Now, let A denote the semi-infinite upper shift 



matrix. Then we observe that multiplying the measure da{x) by x^ or, multiplying d/3(y) by y-', is 
tantamount to multiplying I on the left by A*, or on the right by (A"^)^ respectively, which gives us a 
whole family of bimoment matrices associated with the same K{x, y) but different measures. Thus we 
have 

Corollary 2.1. For any nonnegative integers i,j the matrix of generalized bimoments A*X(A"^)-' is TP. 
We conclude this section with a few comments about the scope of Theorem 12.11 



Remark 2.2. Provided that the negative moments are well defined, the theorem then applies to the doubly 
infinite matrix i,j € Z. 

Remark 2.3. // the intervals are R and K{x,y) = e^^ then the proof above fails because we cannot 
re-define the measures by multiplying by powers of the variables, since they become then signed measures, 
so in general the matrix of bimoments is not totally positive. Nevertheless the proof above shows (with 
a = b = Oora,bCz 2Zj that the matrix of bimoments is positive definite and -in particular- the 
biorthogonal polynomials always exist, which is known and proved in J 15^ . 



2.1 Biorthogonal polynomials 

Due to the total positivity of the matrix of bimoments in our setting, there exist uniquely defined two 
sequences of monic polynomials 



Pnix) = X" + 



such that 



Pn{x)qm{y)K{x,y)da{x)d(3{y) = 
Standard considerations (Cramer's Rule) show that they are provided by the following formulae 

^00 ■ • ■ lon 



Pn{x) 



1 



■ det 



Ion-1 



1 



9n(2/) 



1 



■det 



, Dn+1 ^ „ 
hn = > 0, 



(2-6) 
(2-7) 
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where Dj > by equation (|2-5p . For convenience we re-define the sequence in such a way that they are 
also normalized (instead of monic), by dividing them by the square root of /i„; 



<in{y) 



:det 



00 



Ion-1 



:det 



iOO 



In- 
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In— In 



y 



(2- 



(2-9) 



Thus (p„|g„ 



We note also that the BOPs can be obtained by triangular transformations of [x] , [y] 
p{x)=Sp[x], q(y) = 5,[y], [x] = [l,x,x'^ , . . -Y 



(2-10) 



where Sp.q are (formally) invertible lower triangular matrices such that Sp^{Sq^)'^ = X, where, we 
recall, 2 is the generalized bimomcnt matrix. Moreover, our BOPs satisfy, by construction, the recursion 
relations: 

xp.i{x) = Xi^i+ip^+i{x) + X^^iPi{x) -\ XifiPo{x), yqi{y) = Yi.i+iqi+i{y) Yi^iq^{y) ^ yi,o9o(y), 

which will be abbreviated as 

xp(x) = Xp(a;) , yci{yf = q(2/)^Y^, (2-11) 

where X and Y are Hessenberg matrices with positive entries on the supradiagonal, and p(a;) c[{y) are in- 
finite column vectors p(a;)^ := {pQ{x),pi{x),p2{x), . . .)\ q(y)^ := {qo{y), qi{y), q2{y), ■ ■ - V respectively. 

The biorthogonality can now be written as (p|q"^) = Id where Id denotes the semi-infinite identity 
matrix. Moreover 

(a;p|q^) = X , (pj^q^) = (2-12) 

Remark 2.4. The significance of the last two formulas lies in the fact that the operator of multiplication 
is no longer symmetric with respect to the pairing (•!•) and as a result the matrices X and 'Y'^ are 
distinct. 

2.2 Simplicity of the zeroes 

In this section we will use the concept of a Chebyshev system of order n and a closely related concept of a 
Markov sequence. We refer to |23| and [17j for more information. The following theorem is a convenient 
restatement of Lemma 2 in [T7], p. 137. For easy display we replace determinants with wedge products. 
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Theorem 2.2. Given a system of continuous functions {ui{x)\i — Q ■ ■ ■ n\ let us define the vector field 

\i{x) = [uq{x), u„(a;)]^, x&U. (2-13) 

Then {ui[x)\i = • • • n} is a Chebyshev system of order n on U iff the top exterior power 

u(xo) A u(xi) A • • •u(a;„) 7^ (2-14) 

for all Xq < Xi < ■ ■ ■ < Xn in U . Furthermore, for {ui{x)\i = • • • }, if we denote the truncation of u{x) 
to the first n + 1 components by u„(x), then {ui{x)\i ~ - • •} is a Markov system iff the top exterior 
power 

Un{xo) A u„(a;i) A • • • u„(x„) ^ (2-15) 

for all xo < xi < ■ ■ ■ < Xn in U and all n G N. 

The following well known theorem is now immediate 

Theorem 2.3. Suppose {ui{x)\i = - • - n} is a Chebyshev system of order n on U, and suppose we are 
given n distinct points U. Then, up to a multiplicative factor, the only generalized polynomial 

P{x) = X]r=o '^i'^ii^); which vanishes precisely at Xi, ■ ■ ■ Xn in U is given by 

P(a;) u(a;) A u(a;i) A • • •u(a;„) (2-16) 

Theorem 2.4. Denote by Ui(x) ~ J K(x, y)y^ d(3{y) , i = ■ ■ ■ n. Then {ui(x)\i = ■ ■ - n} is a Chebyshev 
system of order n on M+. Moreover, P{x) as defined in Theorem \2.S\ changes sign each time x passes 
through any of the zeros Xj . 

Proof. It is instructive to look at the computation. Let xo < xi < ■ ■ ■ Xn, then using multi- linearity of 
the exterior product, 

P{xq) = u(a;o) A u(a:i) A • • • u(a;„) = 

K{xo,yo)K{xi,yi) ■ ■ ■ K{x„,y„)[yQ]n A [yi]„ A • • • A [yn]„df3{yo) ■ ■ ■ df3{yn) = 



1 / det[K{x,,yj)]l^^o^iY)df3iyo) ■ • •d/?(y„) = / det[if(x„ %)]^^-=o A(y)d/3(yo) • • •d/?(yn), 

where [y]„ ~ [?/°, y^, ... y"]^ . Thus /'(.tq) > 0. The rest of the proof is the argument about the 
sign of the integrand. To see how sign changes we observe that the sign of P depends only on the ordering 
of x,xi,X2, - ■ ■ Xn, in view of the total positivity of the kernel. In other words, the sign of P is sgn^n) 
where tt is the permutation rearranging x, xi, X2, ■ ■ • Xn in an increasing sequence. □ 

Corollary 2.2. Let {/^(x) -.^ J K{x,y)qr{y)dp{y),\i = O---}. Then {fiix)\i = • • • n} zs a Markov 
sequence on , 
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Proof. Indeed, Theorem l2.2l implies that the group GL{n+l) acts on the set of Chebyshev systems of order 
n. It suffices now to observe that qj are obtained from {1, y, - ■ ■ , y"} by an invertible transformation. □ 

Remark 2.5. Observe that {fi{x)\i = - • - n} is a Markov sequence regardless of biorthogonality. 

Biorthogonahty enters however in the main theorem 

Theorem 2.5. The zeroes ofpn,qn are 0,11 simple and positive. They fall within the convex hull of the 
support of the measure da (for pn 's) and d/3 (for the qn 's). 

Proof. We give first a proof for p„. The theorem is trivial for n = 0. For 1 < ^t- , let us suppose p„ has 
r < n zeros of odd order in the convex full of supp{da). In full analogy with the classical case, 1 < r, 
since 

J Pn{x)fo{x)da{x) = jj Pn{x)K{x,y)da{x)d(3{y) = 

by biorthogonality, forcing, in view of positivity of K{x,y), Pn{x) to change sign in the convex hull 
of supp{da). In the general case, denote the zeros by xi < 2:2 < • ■ -Xr. Using a Chebyshev system 
fi{x), i = 0, ■ ■ ■ r on R+ we can construct a unique, up to a multiplicative constant, generalized polynomial 
which vanishes exactly at those points, namely 

R{x) = F{x) A F{xi) A F{x2) A • • • A F{xr) (2-17) 

where 

F{x)=[fo{x) f,{x) ■■■ fr{x)]\ XGR. 

It follows then directly from biorthogonality that 

Pn{x)F{x) A F{xi) A F{x2) A • • • A F{xr)da{x) = (2-18) 

On the other hand, R{x) is proportional to P{x) in Theorem 12.31 which, by Theorem 12.41 changes sign 
at each of its zeroes, so the product Pn{x)R{x) is nonzero and of fixed sign over R+ \ {xi,X2, ■ ■ ■ ^Xr}. 
Consequently, the integral is nonzero, since a is assumed to have infinitely many points of increase. Thus, 
in view of the contradiction, r > n, hence r = n, for p„ is a polynomial of degree n. The case of qn follows 
by observing that the adjoint K* is also a TP kernel and hence it suffices to switch a with (3 throughout 
the argument given above. □ 

Lemma 2.1. In the notation of Corollary f nix) has n zeros and n sign changes in the convex hull 
of supp{da). 
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Proof. Clearly, since {ui{x)\i ~ ■ ■ - n} is a Chebyshev system of order n on R+, the mmiber of zeros of 
/„ cannot be greater than n. Again, from 

fn{x)po{x)da{x) = 0, 

we conclude that /„ changes sign at least once within the convex hull of supp{da). Let then xi < X2 < 
■ ■ ■ Xr, 1 < r < n he all zeros of /„ within the convex hull of supp(da) at which /„ changes its sign. Thus, 
on one hand, 



/ ^ ~ Xi)fn{x)da{x) > 0, 

i=l 

while, on the other hand, using biorthogonality we get 

j ^ ~ Xi)fn{x)da{x) = 0, 



e = ±, 



which shows that r = n. □ 

In view of Theorem 12.31 the statement about the zeros of /„ has the following corollary 
Corollary 2.3. Heine-like representation for /„ 

fn{x) = Cu{x) A u{xi) A u{x2) • • • A u{xn) (2-19) 
where Xj are the zeros of /„ and C is a constant. 



3 Cauchy BOPs 

From now on we restrict our attention to the particular case of the totally positive kernel, namely, the 
Cauchy kernel 

K{x,y) = ^ (3-1) 
X + y 

whose associated biorthogonal polynomials will be called Cauchy BOPs . Thus, from this point onward, 
we will be studying the general properties of BOPs for the pairing 

//p.(..).™(y)^^^^^^ = (p,J.,„)^ (3-2) 
J J x + y 

Until further notice, we do not assume anything about the relationship between the two measures da, d/3, 
other than what is in the basic setup of Definition 12.51 
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3.1 Rank One Shift Condition 

It follows immediately from equation (j3-ip that 

which, with the help of the shift matrix A and the matrix of bimomcnts X, can be written as: 

a = (ao, Oil, • ■ • 7 '^j - 



(3-3) 



/3 = (/3o,/3i,...)^ , Pi 



x^da{x) > 0, 
y^dl3{y) > 0. 



Moreover, by linearity and equation (j2-12p . we have 

X+Y^ = 7r?7^, TV := J pda , q := J qdp , p(x) := (po(a;),pi (a;), ...)*, q(j/) (^oCy), gi(y), ■•■ )* 

(3-4) 

which connects the multiplication operators in Ha and Hfj. Before we elaborate on the nature of this 
connection we need to clarify one aspect of equation 



Remark 3.1. One needs to exercise a great deal of caution using the matrix relation given by equation 
P-4p . Its only rigorous meaning is in action on vectors with finitely many nonzero entries or, equivalently, 
this equation holds for all principal truncations. 

Proposition 3.1. The vectors TT,ri are strictly positive (have nonvanishing positive coefficients). 

Proof. We prove the assertion only for tt, the one for t] being obtained by interchanging the roles of da 
and d/3. 

From the expressions (|2-9I) for Pn{x) we immediately have 



1 



DnDn+l 



■ det 



^00 



ao 



(3-5) 



Since we know that Z3„ > for any n > we need to prove the positivity of the other determinant. 
Determinants of this type were studied in Lemma 4.10 in [21| . 

We nevertheless give a complete proof of positivity. First, we observe that 

. n+l 



/ d"+ a 
n^?''n^r' n" (r 



A(xri)n; 



_i d"+iad"/3 
]Ti=i{xj +yj)' 



(3-6) 
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Here the symbol X"~^^ is to remind that the vector consists of n + 1 entries (whereas Y consists of n 
entries) and that the Vandermonde determinant is taken accordingly. Note also that the variable Xn+i 
never appears in the product in the denominator. Symmetrizing the integral in the Xj's with respect to 
labels j = 1, . . . , n , but leaving Xn+i fixed, gives 

d"+iad"/3 



J n,=i( 



Symmetrizing now with respect to the whole set xi, . . . , Xn+i we obtain 

K{xi,yi) ... K{xn+i,yi) 

1 ' 



(3-7) 



n!(n + 1)! 



A{Xl'+')A{Y)dct 



K{xi,yn) ... K{Xn+l,yn) 

1 ... 1 



d"+^ad"/3 (3-8) 



Moreover, since the integrand is permutation invariant, it suffices to integrate over the region < 
xi < X2 < ■ ■ ■ < Xn < Xn+i X < i/i < 2/2 < ' ' ' < J/n, and, as a result 



TTn\/Dn+lDn 



0<xi<X2<-<Xn + iO<yi<y2<-<yn 



K{xi,yi) ... K{xn+i,yi) 

K{xi,yn) ... K{Xn+l,yn) 

1 ... 1 



d"+iad"/3. 



(3-9) 



We thus need to prove that the determinant containing the Cauchy kernel is positive for < xi < 
X2 < ■ ■ ■ < Xn+i and < yi < y2 < ■ ■ ■ < yn- It is not difficult to prove that 

1 1 



det 



xi+yi 



xi+yn 
1 



Xn+i+yi 



Xn+l+ym 

1 



A(Xi"+^)A(y) 

Ki^lUl^ii^j +yk) 



(3-10) 



and this function is clearly positive in the above range. 



□ 



3.2 Interlacing properties of the zeroes 

From ([2]), (|2-10|1 and (|2-lip the following factorizations arc valid for all principal truncations: 

Moreover, since T is TP, the triangular matrices and S'-^ are totally nonnegative (TN) [13] and have 
the same diagonal entries: the nth diagonal entry being D„/D„_i. Furthermore, one can amplify the 
statement about S'^ and using another result of Cryer (|12j) which implies that both triangular 
matrices are in fact triangular TP matrices (all non-trivial in the sense defined in [12] minors are strictly 
positive). This has the immediate consequence 



15 



Lemma 3.1. All principal truncations 'K.[n],Y[n\ are invertible. 



Proof. From the factorization X = SpA{Sp) ^ we conclude that it sufhces to prove the claim for AS^ ^ [n] 
which in matrix form reads: 



iSp^)20 iSp^)21 iSp^)2 



{Sp 



i^p ^)n+l,0 (Sp 



+ 1,1 







+ l,n 



□ 



However, the determinant of this matrix is strictly positive, because S'p ^ is a triangular TP. 

Remark 3.2. This lemma is not automatic, since A[n] is not invertible. 

We now state the main theorem of this section. 
Theorem 3.1. X and Y are TN. 



Proof. We need to prove the theorem for every principal truncation. Let 72 > be fixed. Wc will suppress 
the dependence on n, for example X in the body of the proof means X[n] etc. First, we claim that X 
and Y admit the L-U factorization: X = X_X+, Y = Y_Y-|_, where denotes the upper triangular 
factor and is the unipotent lower triangular factor in the Gauss factorization of a matrix A. Indeed, 
X_|_ = (AS'~^)+, Y_|_ = (A5'~^)+ are upper triangular components of TN matrices AS~^ and AS~^ and 
thus are totally nonnegative invertible bi-diagonal matrices by Lemma |3. II 
From X + Y-^ — ttt]^ we then obtain 



(YT)-iX_ + Y_X;i = {{Ylr^n) (ry^X-i^ 



We need to show that vectors p , p, have positive entries. For this, notice that 

p = {{Y+fr'Spc = {{{As-'urr'Spc , 

p - {{X+f)-'S,f3 = {{{As;')+f)-'s,i3. 

Now, it is easy to check that if the matrix of generalized bimoments I is replaced by TA^ (see 
Corollarv l2.ll ') then Sp {{{AS~^)+)'^)~^ Sp, while a is unchanged, which implies that p is a new tt in 
the notation of Proposition 13.11 and hence positive by the same Proposition. Likewise, considering the 
matrix of generahzed bimoments AX, for which f3 is unchanged, Sq {{{ASp^)+)'^)^^Sq and /j is a new 
r] in the notation of Proposition 13.11 implying the claim. 

Thus 

p = Dpi ,p = Dp^l, 
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where Dp , are diagonal matrices with positive entries and 1 is a vector of Is. 
We have 

D-\Yl)-'^^D-' + D-^Y^y.-'D-' = 11^. 

The first (resp. second) term on the left that we can call X (resp. Y"^) is a lower (resp. upper) triangular 
matrix with positive diagonal entries . The equality above then implies that (i) Xij = Yij = 1 for all 
i > j and (ii) Xa + Ya ~ 1 for all i. In particular, both Xa and Ya are positive numbers strictly less 
then 1. 

This means that X, Y admits factorizations 

±={Id- K^y^Lx , Y = {Id- K^)-^Ly , 

where 

oo oo 

Lx = ^ XiiEii + (1 — Xii)Ei+i i , Ly — ^ YiiEii + (1 — Yii)Eij^i i . 

1=0 i=0 

Since all entries of bi-diagonal matrices Lx, Ly are positive, these matrices are totally nonnegative and 
so are 

X = Yl{Id - A^)-iLxX+ , Y = X^(/d - K^)-'^LyY+ . (3-11) 

□ 

Corollary 3.1. X and Y are oscillatory matrices. 

Proof. We give a proof for X. The factorization (|3-11[) we have just obtained shows that X is the product 
of an invertible lower-triangular TN matrix Y!^(/rf — A-^)^^ and a tri-diagonal matrix J = La'X+. Note 
that Lx has all positive values on the main diagonal and the first sub-diagonal. Entries on the first super- 
diagonal of X+ coincide with corresponding entries of X and thus are strictly positive by construction. 
Moreover, leading principal minors of X arc strictly positive (see the proof of Lemma l3.1|) . which implies 
that all diagonal entries of X+ arc strictly positive too. Thus J is a tri-diagonal matrix with all non-trivial 
entries strictly positive. 

Since diagonal entries of Y!^(/(i — A-^)^^ are strictly positive and all other entries are non- negative, 
every zero entry of X implies that the corresponding entry of J is zero. In view of that all entries on 
the first super- and sub-diagonals of X must be strictly positive, which, by a fundamental criterion of 
Gantmachcr and Krcin (Theorem 10, II, [13), ensures that X is oscillatory. 

□ 

The interlacing properties for the zeros of polynomials p„,g„, as well as other properties of Sturm 
sequences, follow then from Gantmacher-Krein theorems on spectral properties of oscillatory matrices (see 
II, Theorem 13, in [17j). We summarize the most important properties implied by Gantmacher-Krein 
theory. 
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Theorem 3.2. The sequences of BOPs {qn} and {pn} are Sturm sequences. Moreover, 

1. their respective zeros are positive and simple, 

2. the roots of adjacent polynomials in the sequences are interlaced, 

3. the following alternative representations of the biorthogonal polynomials hold 



Pn{x) ^ \ J^'^ det{x ~ X[n- 1]), 1 < n, 

V -t^n+l 



Remark 3.3. The fact that the roots are positive and simple follows indeed from the fact that X and Y 
are oscillatory. Theorem p.Sp . however, indicates that this property is true even for a more general case 
when the totally positive kernel K(x,y) is not necessarily the Cauchy kernel. 

4 Four-term recurrence relations and ChristofFel Darboux iden- 
tities 

We establish in this section a basic form of recurrence relations and an analog of classical Christoffel- 
Darboux identities satisfied by {qn} and {pn}- First, we introduce the following notation for semi-infinite, 
finite-band matrices. 

Definition 4.1. Given two integers a <b , a semi-infinite matrix A is said to have the support in [a,b] 

j — i < a or j — i > b imply Aij = (4-1) 
The set of all matrices with supports in [a, b] is denoted M[a,b] ■ 

The content of this section relies heavily on the relation (|3-4[) which we recall for convenience: 



where Dt^, D,, respectively, are diagonal matrices of averages of p and q. Since the vector 1 is a null 
vector of A — /d we obtain 

Proposition 4.1. X and Y satisfy: 

1. (A - Id)D-^X + (A - Id)D-^Y'^ = 0. 

2. A:= (A-/d)i5-iXeM[_i^2]- 
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3. Xi?-i(A^ - Id) + Y^D-i(A^ - Id) = 0. 

4. A:^y.D-\KT -Id) e Af[_2,i]. 

As an immediate corollary we obtain the factorization property for X and Y. 
Corollary 4.1. Let A, A and 

L -.^ {A- Id)D-\ L := L»^i(A^ -/d), 
respectively, denote matrices occurring in Proposition \4- 1\ Then 

LX^A, xl^A, ^ e A/[_i.2], A e Af[_2,i]- 
Likewise, Y admits a similar factorization: 

YL^ = B, (l'^)Y = B, 

where B = -A^, B = -A^ . 
Hence, 

Corollary 4.2. p and q satisfy four-term recurrence relations of the form 

'Pn{x) Pn~i{xy 



lPn+l{x) + An-l^nPn{x) + A„_i,„_ip„_i (x) + A„_i,„_2Pn-2 (-t) , 



y I Qjiiy) _ ^" '^^y^ \ ^ S„_i_„+iq„+i(j/) + Bn-i,„qniy) + S„_i_„_ig„_i(y) + Bn-l,n-2qn-2{y), 



Vn Vn-l 

for I < n with the proviso that p_i — q-i =0. 
Proof. Wc give the proof for p(a;) in matrix form. Indeed, from 

xp{x) ~ Xp(a;), 

it follows that 

xLp{x) = LXp(a;), 

hence the claim, since L £ M[o.i] and LX = A E 7\/[_i_2]- □ 

Let us observe that L has a unique formal inverse, represented by a lower triangular matrix. Let us 
then define 

p(a;) = L-^p{x). 
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Theorem 4.1 (ChristofFcl-Darboux Identities for q and p). 

n-l 

{x + y)Y^ qjiy)pj{^) - q^(j/)[n, {y - y^)l]p{x) (4-2) 

where 11 := n„ is the diagonal matrix diag{l, 1, . . . , 1, 0, . . . ) with n ones (the entries are labeled from 
to n — 1). The explicit form of the commutators is: 

[n, (y — Y )L] = yl„_i.„£'„„i.„ — ( h ^„.„_i)iJ„.„_i — 

An,n-2En,n~2 ^ ^ri+l,?i-l-£'n+l,?i-l i (4-3) 

where Ai_j, Ai_j respectively, denote the {i,j)th entries of A, A, occurring in Proposition\4. 1\ 



Proof. We give the proof of equation (|4-2p . Since {y — Y)q = it suffices to prove that the left hand side 
equals q^n(y — Y'^)Lp{x). From the definition of p and equation (|2-Iip we obtain 

{x + 2/)q^(y)np(.T) = q^(2/)nyip(x) + q^(2/)nXp(x) = q^(2/)n2/Lp(x) + q^(y)nXLp(x), 

which, after switching XL with — Y^L in view of Proposition 14.11 gives equation (|4-2p . To get the 
commutator equation (|4-3[) one needs to perform an elementary computation using the definition of 
A. □ 

We establish now basic properties of p and its biorthogonal partner q defined below. 

Proposition 4.2. The seguences of polynomials 

p = L-ip , q^ = q^L (4-4) 

are characterized by the following properties 

1. degg„ degp„ = n; 

2. I g„d/3 = 0; 



3- II Pn(a:)gm(y) ^"^'^ = ^mn / 

a; -I- ?/ 



I g„(2/) = ;^y^^y"+i+0(2/"); 
/n addition 

a. q and p satisfy the intertwining relations with q and p 

yq^ - -q'^A 

xp = Ap; (4-5) 
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b. q and p admit the determinantal representations: 



9n(y) 



VnVn+l \jDnDn+2 



: det 



Pn{,x) 



■ dct 



^00 



InO 



-'00 

In-10 
1 

-^n— 1 Ti 
n 



^ 0)1+1 



yU+l 



(4-6) 







(4-7) 



^ , , dad/3 
x + y 



Proof. Assertions (1), (2) and (4) follow directly from the shape of the matrix L. Assertion (3) follows 
from (p,q') = 1 by multiplying it by L on the right and by L^^ on the left. Assertion (c) follows from 
assertions (1), (2) and (3); indeed from (2) and (3), it follows that the polynomial p„ is biorthogonal 
to all polynomials of degree < n with zero d/3-avcrage and {PoV'' ~ Pj ^ < j < n} is a basis for such 
polynomials. 

The intertwining relations follow from the definitions of the matrices L, A and of the polynomials p, q. 
The determinantal expression for q„ follows by inspection since the proposed expression has the 



defining properties (1) and (2) and is biorthogonal to all powers l,x, 



The normalization is 



found by comparing the leading coefficients of g„ = -^^^^qn+i + 0{y"). The determinantal expression for 



Pn{x) follows again by inspection; indeed if F{x) is the determinant in (j4-7p then 

^00 ■ • • ^0 ri loj 



{F{x)W)=Act 



InO 

Po 



—In — 1 j 

In 71 j 
Pn 



= ~P,D^+^^^{F{x)\l 
Po 



(4-8) 



where the determinants are computed by expansion along the last row. The proportionality constant is 
again found by comparison. □ 



One easily establishes a counterpart to Theorem 14. II valid for q and p. 
Proposition 4.3 (Christoffcl-Darboux identities for q and p ). We have 



{x + y)Y, %iy)pjix) - q^(2/)[(a; - X)£,n]p(x) - q^(y)[n, {-x - F^)L]p(x). (4-9) 
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Remark 4.1. Observe that the commutators occurring in both theorems have identical structure; they only 
differ in the variable y in Theorem \4-.l\ being now replaced by —x. We will denote by A{x) the commutator 
[n, {—X — Y"^)L] and by A„(x) its nontrivial 3x3 block. Thus the nontrivial block in Proposition^ 
reads: 








^n— l.n 




An{x) = 







(4-10) 




—An+l^n-l 








while the block appearing in Theorem \4.1\ is simply A„(— y). 

With this notation in place we can present the ChristofFel-Darboux identities in a unified way. 

Corollary 4.3 (Christoffel-Darboux identities for q, p, and q, p ). The biorthogonal polynomials q, p, 
and q, p satisfy 

n-l 

{x + y)Y. iMpA^) - q^(y)A(-y)p(x), (4-11) 



3=0 
n-l 



{x + y)Y^ %{y)pj{x) = q^{y)A{x)p{ 



(4-12) 



3=0 



5 Approximation problems and perfect duality 

We wiU associate a chain of Markov functions associated with measures da and d/3 by taking the Stieltjes' 
transforms of the corresponding measures as weU as their reflected, with respect to the origin, images. 



Definition 5.1. Define 

1 



V 



-my). 



1 



-da(x)d/3(2/), 



W^. [z) = 
Wpo.' {z) = 



1 



Z + X 



Aa{x) 



1 



-da{x)dp{y). (5-1) 



{z + x){x + y) JJ {z-y){y + x) 

We recall now an important notion of a Nikishin system associated with two measures (see [23], p. 
142, called there a MT system of order 2). 

Definition 5.2. Given two measures d^i and d/i2 with disjoint supports Ai, A2 respectively, a Nikishin 
system of order 2 is a pair of functions 

dAil(a;i) J. f dfii{xi) f d^i2{x2) 

/ai z- 



h{z) 



Hz) 



/Ai Z~Xi Jai Z-Xi J^^Xi-X2 

Remark 5.1. The definition of a Nikishin system depends on the order in which one "folds" measures. 
If one starts from d/i2 , rather than d/ii one obtains a priory a different system. As we show below 
the relation between these two Nikishin systems is in fact of central importance to the theory we are 
developing. 
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The following elementary observation provides the proper framework for our discussion. 

Lemma 5.1. Let da* denote the measure obtained from da by reflecting the support of da with respect 
to the origin. Then Wp, Wpa* (^nd Wa* , W^tp are Nikishin systems associated with measures dfi and da* 
with no predetermined ordering of measures. 

The relation between these two Nikishin systems can now be readily obtained. 

Lemma 5.2. 

Wp{z)W^,{z) ^Wp^.{z) + W^.fi{z). (5-2) 
Proof. Elementary computation gives: 

which implies the claim. □ 

Remark 5.2. Equation (|5-2p was introduced in \21^ for the DP peakons (see Lemma Jf..l there). Observe 
that this formula is valid for any Nikishin system of order 2. 

We formulate now the main approximation problem, modeled after that of |21| 

Definition 5.3. Let n > 1. Given two Nikishin systems Wp^Wpa- and W^'^Wa'p we seek polyno- 
mials Q{z), degQ ~ n, Pp{z), degPp = n — 1 and Ppa*{z),degPj3a-' = n — 1, which satisfy Pade-like 



approximation conditions as z ^ oo, z G C±.' 

Q{z)Wp{z) - Pp{z) = o(^^y (5-3a) 
Q{z)Wo.,p{z) - Pp{z)Wo.,{z) + Pp^.{z) = O (5-3c) 



Remark 5.3. In the case that both measures have compact support we can remove the condition that 
z g C± since all the functions involved are then holomorphic around z = oo. 

Remark 5.4. In the terminology used for example in J2^ the triplets of polynomials Q, Pp, Ppa' provide a 
Hermite-Pade approximation of type I to the Nikishin system Wp^Wfja- and, simultaneously, a Hermite- 
Pade approximation of type II to the Nikishin system Wa* , Wa'p. 

Definition 5.4. We call the left hand sides of approximation problems (|5-3p Rp^Rp (J* and Ra* p respec- 
tively, referring to them as remainders. 



z — y z + X 



da{x)d(3{y), 
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The relation of the approximation problem (|5-3p to the theory of biorthogonal polynomials q and p 
is the subject of the next theorem. 

Theorem 5.1. Let qn{y) be defined as in (|2-9p . and let us set Q{z) = Qniz) Then Q{z) is the unique, 
up to a multiplicative constant, solution of the approximation problem (j5-3l) . Moreover, Pf3,Pj3a' o,nd all 
the remainders Rp, Rpa* and Ra*i3 are uniquely determined from Q with the help of the formulas: 

. / m^m^p^y)^ P,„.(,) = // ^^^l-9M_Mx)my\ (5-4a) 
J z-y J J (z-y)(x + y) 

Rp{z)= [ ^dPiy), Rpo.'{z) = II 7 ^[f ^ M^Wiv), (5-4b) 

J z-y J J [z-y){x + y) 

Rc'piz) ^-JJ -^^^^^^—^AaixWiy) = j ?±^da*{x). (5-4c) 

Proof. We start with the first approximation problem involving Q{z)Wp{z). Writing explicitly its first 
term we get: 

OW.a,,fa,./0<£)^d«,,+ /iMa«,,. 



Z-y J z~y J z-y 

Since J ^^^jz^^dp{y) is a polynomial in z of degree n ~ 1, while J ^^dP{y) — 0{^), we get the 
first and the third formulas. The second and fourth formulas are obtained in an analogous way from the 
second approximation problem. Furthermore, to get the last formula we compute P/j and Ppa* from the 
first two approximation problems and substitute into the third approximation problem, using on the way 
Lemma [521 to obtain: 

RfsWa' — R(3a' = Ra*l3- 

Substituting explicit formulas for Rp and Rpa* gives the final formula. To see that Q{z) is proportional 
to qn{z) we rewrite —Ra*p as: 



-da{x)dP{y)- ^^y^ 



{z + x){x + y) J J (x + y) 



1 i-(-(f)r 



Z + X Z + X 



da{x)d(3{y)- 



V ^'!i\ , dadp^ 11-^, ^ da{x)dp{y)+ // V ^—^^^^ dad/3 



j^o ""^^^ ix + y){z + x) J J (x + y) 



Z + X 



n-1 



{~xy Q{y) 



^-^0 2^+1 {x + y){z + x) 



To finish the argument we observe that the first term is already 0{—^), hence the second term must 
vanish. This gives: 

^^^da(a;)d/3(y) =0, Q<j<n-l, 
x + y 

which characterizes uniquely (up to a multiplicative constant) the polynomial qn- □ 

Remark 5.5. In the body of the proof we used an equivalent form of the third approximation condition, 
namely 

RpWa^ (z) - Rp^, (z) = Ro^-^niz) = O(^). (5-5) 
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By symmetry, we can consider the Nikishin systems associated with measures a and (3* with the 
corresponding Markov fimctions Wa,Wap* and Wp- ,Wi3*a- We then have an obvious interpretation of 
the polynomials p„ . 

Theorem 5.2. Let pn{x) be defined as in (|2-9p . and let us set Q{z) ~ pn{z). Then Q{z) is the unique, 
up to a multiplicative constant, solution of the approximation problem for z oo, z G <C±: 

g(z)VP^„(z)-F„(z) = , (5-6a) 

Q{z)W^f,. (z) - P^f,. (z) = <^ (^) ' (5-6b) 

Q{z)Wp,a.iz)-Pc.{z)Wp,{z) + P^p,{z) = o(^-^^ , (5-6c) 

where Pa, Pa/3* are given by formulas of Theorem \5.1\ after switching a with [3. 

Clearly, one does not need to go to four different types of Nikishin systems in order to characterize 
qn and p„ . The following corollary is an alternative characterization of biorthogonal polynomials which 
uses only the first pair of Nikishin systems. 

Corollary 5.1. Consider the Nikishin systems Wp, Wpa* and Wa* , Wa'p- Then the pair of biorthogonal 
polynomials {g„,p„} solves: 

1. Q{z) = qn{z) solves Hermite-Pade approximations given by equations (|5-3p . 

Q{z)Wp{z)-Pf3{z)=o(- 



Qiz)Wi3a'{z)~Ppa'iz)=0 
Qiz)Wa'piz) - P[3iz)Wa^{z) + Ppa'{z) - O 



2. Q{z) ~ pn{—z) solves switched (Type I with Type II) Hermite-Pade approximations 

Q{z)Wa'{z)-Pa*{z)=o(-] , (5-8a) 



Q{z)Wa^l3{z) - Pa'p{z) = O I - ) , (5-8b) 



z 

Q{z)W0a'{z) - Pa'{z)Wp{z) + P„.^(z) = O (^^) (5-8c) 

We finish this section with a few results needed for the Riemann-Hilbert problem approach to biorthog- 
onal polynomials {<7„, which will be presented in the next section. 
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Definition 5.5. We define the auxiliary vectors in addition to the main polynomial vectors cig{vu) := c{{w) 
and P(, (z) := p{z), as 



iM := J niy) 



q 



my) 

w-y' 
p{x)da{x) 

Z — X 



w — X 

v>i{y) 

z-y 



da*(x), 



d/3*(2/). 



Moreover, 



Here 1 is the vector of ones 



P2(^) := j 



1 



p,(z) := p,(z) + -(p|l) = i-ip,(z) - 1, 



Pi(?/) 



d/3*(y). 



(5-9) 
(5-10) 

(5-11) 
(5-12) 



Remark 5.6. Note that the definition above unifies the approximants and their respective remainders 
(see Theorem \5.1] ). thus, for example, qi(u') = R^(u'), (w) ~ 'R,a*i3{w) etc. The definition of "hatted" 
quantities is justified below. 



Theorem 5.3 (Extended Christoffel-Darboux Identities). Let a, 5 = 0, ... 2. Then 

{w + z)q^(w)np^ (z) = q^(w)A(-u;)p^ (z) - ¥{w, z)ab 

where 

¥{w,z) 



(5-13) 



1 



(5-14) 





1 Wp*{z) + Wp{w) 

1 Wo.{z) + Wo.'{w) Wo.'{w)Wp,{z) + Wc.*p{w)+Wp*c.{z)_ 

Proof. The proof goes by repeated applications of the Christoffel-Darboux Identities given by Theorem 
14.11 and Pade approximation conditions 15-31 The details have been relegated to Appendix El □ 



We point out that if we set w — — z in the GDI's contained in Theorem 15.31 the left hand side 
vanishes identically and the RHS contains terms of the form qa(— z)A(z)pj,(z) minus Fq(,(— z,z). The 
main observation is that the second term is constant, independent of both z and n, and hence one ends 
up with the perfect pairing (see [3]) between the auxiliary vectors. For the reader's convenience we 
recall the definition of A(z) to emphasize the implicit dependence on the index n hidden in the projection 

n. 



Theorem 5.4. (Perfect Duality) 
Let 



1 

1 

1 



Then 



qr(- 



(^)Pb(^) = Jafc, 



where k{z) = [(z - X)L,n]. 



'The formula /3q ' < p„, 1 >= —1 follows directly from the determinantal expression in Proposition [42] 
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Proof. The only nontrivial entry to check is (2,2). In this case, after one substitutes w = ~z into 
Wa_' {w)Wf}' (z) + Wa*si{vS) + W/3'a{z), One obtains the identity of Lemma [5^ □ 



There also exists an analog of the extended Christoffel-Darboux identities of Theorem 15.31 for the 
"hatted" quantities. 
We first define: 

Definition 5.6. For a = 0, 1, 2, 

:= q^L. (5-15) 
The following identities follow directly from the respective definitions. 
Lemma 5.3. 

\liw)Yn, a = 0,1 

q^{w)Yn-{l\iS), a = 2. 

b = 0, 
5=1, 

Po|l) + ^2^^±f^, b = 2. 
Theorem 5.5 (Extended Christoffel-Darboux Identities for qa, Pb). Let a, 6 = 0, ... 2. Then 

{w + z)ql{w)Upb{z) = ql{w)A{z)pb{z) - ^{w, z)ab (5-16) 

where 




V{w,z) =¥{w,z) ""^^ 



1 W/5.(z) 

Wp{z) Wpiw)Wp,{z) 

1 Wo^.p.iw) Wo,.0.{w)W0.{z) 



Proof. We give an outline of the proof. For a = 0, 1, in view of Lemma [ 

(w + z)q^(w)npf,(z) = q^(w)A(z)pt(z) + q^(w)n(z - X)Lpb(z). 
The second term equals, again by Lemma |5.3[ 



(5-17) 



0, 6 = 0, 

\z + 



qr(^)n^^2^, 5=1, 



Now, one goes case by case, using biorthogonality of q^ and po, and the definition of q[{w). After a 
few elementary steps one arrives at the claimed result. The computation for a = 2 is only slightly more 
involved. From Lemma 15.31 we obtain: 



{w + z)q^{w)npb{z) = q^{w)A{z)pb{z) - (l|qo)npb(z) + q^ {w)U{z - X)Lpfc(z). 
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In view of biorthogonality of and p, after some intermediate computations, one obtains: 



1, 6 = 



(l|qo)np6(z) 




Likewise, 



qjHn(z-X)Lpb(z) 








1, 

2, 



and the claim follows. 



□ 



6 Riemann— Hilbert problems 

In this section we set up two Riemann-Hilbert problems characterizing the Cauchy BOPs that enter the 
Christoffel-Darboux identities of the previous section. This is done in anticipation of possible applications 
to the study of universality for the corresponding two-matrix model. Moreover, since the Christoffel- 
Darboux kernels contain also the hatted polynomials, it is useful to formulate the Riemann-Hilbert 
problems for those polynomials as well. 

We will also make the assumption (confined to this section) that the measures da, d/3 are ab- 
solutely continuous with respect to Lebesgue's measure on the respective axes. Thus one can write 
^ = e s"^ , 1^ = e for the respective (positive!) densities on the respective supports: the 

signs in the exponents are conventional so as to have (in the case of an unbounded support) the potentials 
U, V bounded from below. The constant h is only for convenience when studying the asymptotics of 
biorthogonal polynomials for large degrees (small h) . 

Since the Christoffel-Darboux identities involve the expressions q^ Ap^ , we are naturally led to char- 
acterize the sequences q and p. However, the other sequences can be characterized in a similar manner 
by swapping the roles of the relevant measures and symbols. 

6.1 Riemann-Hilbert problem for the q— BOPs 

We will be describing here only the RHP characterizing the polynomials qn {y) , where the characterization 
of the polynomials Pn {x) is obtained by simply interchanging a with p (see for example Theorem 15. 2p . 
We consider the real axis M oriented as usual and define 
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For simplicity of notation wc will suppress the superscript in most of the following discussions, only 
to restore it when necessary for clarity; the main point is that an arrow on top of the corresponding 
vector will denote a "window" of three consecutive entries of either the ordinary vector q (index a = 0), 
or the auxiliary vectors (index a = 1, 2, see Def. 15. 5p which, as we might recall at this point, combine 
the polynomials and the corresponding remainders in the Hermite-Pade approximation problem given by 
Theorem 15. II Some simple observations are in order. The vector q^ (w) is an analytic vector which has a 
jump-discontinuity on the support of d/3 contained in the positive real axis. As w — * oo (away from the 
support of d/3) it decays as ^. Its jump-discontinuity is (using Plemelj formula) 



it u 



qi(ti;)_ — 27ri^qo(w) , w € supp{dP). 



['nn-2 Vn-l 77„]* + ^(l/w^) 



(6-2) 
(6-3) 



Looking at the leading term at w = cxd we see that 

1 

w 

The vector q^ {w) is also analytic with a jump discontinuity on the reflected support of da (i.e. on 
supp(da*)). In view of Theorem (STTJ recalling that q2 are remainders of the Hermite-Pade approximation 
problem of type II, we easily see that 

t 



C„-l 



The jump-discontinuity of is 



[l + Oil/w)), 



.da* 



D 



n+l 



> 0. 



(i^iw)^ = q2(u')_ — 27ri--^^ — Qi (^) ^ G supp{da*). 



The behavior of q^ {w) at infinity is 

with the same c„'s as in 
Define the matrix 



C„-2 



Cn-l 



W 



(i + o(iA«)), 



(6-4) 



(6-5) 



(6-6) 



1 


(-1)'^ 



1 






(-1)" 












q(")(w;),qW(^),cjW(, 



(6-7) 



Proposition 6.1. The matrix T{'w) is analytic onC \ {supp{dp) U supp{da*). Moreover, it satisfies the 
jump conditions 



1 



1 
1 




-27ri^ 



-27ri 




1 

da* 



w e supp{dl3) C R 
w e supp{da*) C 



(6-8) 
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and its asymptotic behavior at w ^ oo is 

T{w) = [l + Oiw-')) 
Moreover, T{u}) can be written as: 





w-'^ 




(6-9) 











1 


















<?«-! qi,n-l 92,71-1 
qn-1 qi,n-l 92,71-1 
qn-2 qi,n-2 92,71-2 



Proof. All the properties listed are obtained from elementary matrix computations. 



(6-10) 



□ 



Remark 6.1. An analogous problem with the roles of a,f3, etc., interchanged, characterizes the monic 
orthogonal polynomials pn-i{x) of degree n — 1 in x. 

Corollary 6.1. Given n G N, the absolutely continuous measures d/3 C K+ and da* C and assuming 
the existence of all the bimoments lij there exists a unique matrix T{w) solving the RHP specified by 
equations (|6-8p . (j6-9p . The solution characterizes uniquely the polynomials Qn-i as well as Qn-i- In 
particular, the normalization constants c„_i,77„_i (i.e. the "norm" of the monic orthogonal polynomials 
and the /3 average of the qn-i) are read off the following expansions 



r2,i(w) = 



Cn-lVn-l 



or, equivalently. 



1 



Vn-l 



(-1)" lim wT2,iH^2,3H, 



-1)" lim w 



0{w- 



2n-ir2,3 



H 



(6-11) 



(6-12) 



r2,i(w) 

Proof. Given d/3 and da* it suffices to construct the Nikishin systems Wp, Wpa* and Wa' , Wa-'fs followed 
by solving the Hermite-Padc approximation problems given by equations (j5-3p . The existence of the 
solution is ensured by the existence of all bimoments lij (see equation (|2-4p for the definition). Then 



one constructs the polynomials qj, finally the matrix T{w) using equation (j6-10p . By construction T{w) 
satisfies the Riemann-Hilbert factorization problem specified by equations (|6-8p and (|6-9p . Since the 
determinant of T{w) is constant in w (and equal to one), the solution of the Riemann-Hilbert problem 
is unique. The formulas for 7y„_i and c„_i follow by elementary matrix computations. □ 



Remark 6.2. By multiplication on the right with a diagonal matrix Y{w) := r(u')diag ^ exp ^- 



2V+U* 

3h 



exp 



2U*+V 



-g^r — , . one can reduce the RHP to an equivalent one with constant jumps. R then 

follows that Y{w) solves a linear ODE with the same singularities as V',U*'; for example if U' ,V' are 
rational functions then so is the coefficient matrix of the ODE and the orders of poles do not exceed 
those ofV', U' . In this case it can be shown that the principal minors of the matrix of bimoments are 
isomonodromic tau-functions in the sense of Jimbo-Miwa-Ueno USf . 
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6.2 Riemann— Hilbert problem for the p— BOPs 



Referring to the defining properties of pn{x) as indicated in Prop. 14.21 we are going to define a second 
3x3 local RHP that characterizes them. 
Define 

Po(^) := [ Pn-2{z) Pn-liz) p„(z) ]* (6-13) 

and Pi 2{z) as the same windows of the auxiliary vectors pi_2 introduced in Definition 15.51 We first 
study the large z asymptotic behavior of po,7^(2:),Pl,7^(•^),i'2,n(^:)• 



Lemma 6.1. The asymptotic behavior at z —* oo, z G C± is given by: 

PoAz)^- — z"{l + 0{l/z)), 

Pl,„(2) = -1 + 0(1/Z), 



(6-14) 
(6-15) 
(6-16) 



Proof. Wc give a proof for pi^„(z) — J £^^^da(a;) -t- -^{po^n\^). The first term is 0{^), while the second 
term can be computed using biorthogonality and the fact that |3o,n = ^(^MPo.n + fyn-iPo.n-i + • • ■ + '/oJ'o.o)- 
Thus the second term equals ~^(Po,o|l) = ^1, since rjQ = (7o/5o, hence the claim for pi_„(z) follows. The 
remaining statements are proved in a similar manner. □ 

For reasons of normalization, and in full analogy with equation (j6-7p . we arrange the window of all 
ps wave vectors into the matrix 

=■■^fs 



T{z) 













" 1 


-1 


' 





-1 










1 





(-1)" 













-1 


1 













P(^),Pl(^),P2(^) 



(6-17) 



Proposition 6.2. The matrix T{z) is analytic m C \ supp{da) U supp(d(3*). Moreover, it satisfies the 
jump conditions 



r(z)+==r(z)_ 



-27ri^ 

dz 



r(z)+ = r(z)_ 

and its asymptotic behavior at z — oo is 

f(z): 



1 
1 




1 -2TTi^ 
dz 



z e suppiAoL) C M_|_ 
, 2 e supp(d/3*) C R_, 





(6-18) 



1 

^7 



(6-19) 
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r(2:) can he written as: 



r(z) = 



'Cn 

-1 







(-1)- 



PO,n Pl,n P2,n 

P0,n-1 Pl,7i-1 P2ji~l 
PO,n-l Pl,n-1 P2,n-1 



(6-20) 



The existence and uniqueness of the solution of the Riemann-Hilbert problem (|6-18p , (|6-19p is proved 
in a similar way to the proof of Corollary 16. II 

Corollary 6.2. Given n G N, the absolutely continuous measures da C R+ and d/3* C M_, and assuming 
the existence of all the bimoments Iij there exists a unique matrix T{z) solving the RHP specified by 
equations (|6-18p . (|6-19p . The solution characterizes uniquely the polynomials pn^i andpn- 



7 Acknowledgments 

M.B. would like to thank the Department of Mathematics of the University of Notre Dame for hospi- 
tality during which the project was initiated and J. Harnad for insight on the relationship of Cauchy 
biorthogonal polynomials with matrix models. 

While working on this project, M. B. and M. G. enjoyed the hospitality of the Department of Mathe- 
matics, University of Saskatchewan and M. G. and J. S. enjoyed the hospitality of the Centre de recherches 
mathematiques, Universite de Montreal. 

J.S would also like to thank H. Lundmark for an ongoing collaboration on the cubic string problem 
which motivated many of the questions addressed in this paper. 

A Appendix: Proof of Extended Christoffel-Darboux Identities 

Theorem A.l (Extended Christoffel-Darboux Identities 15. 3|) . Let a, 6 = 0, ... 2. Then 

{w + z)qy(w)np, (z) = qJ{w)A{~w)p, (z) - ¥{w, z)ab (1-1) 



where 



¥{w,z) 



1 

1 Wp*{z) + Wi3iw) 

1 WUz) + Wa..iw) W,,,iw)Wf3,iz) + Wa,,p{w)+Wf3.a.iz) 



(1-2) 



Proof. The proof goes by repeated applications of the Christoffel-Darboux Identities given by Theorem 
4.11 and Fade approximation conditions (|5-3p . We observe that all quantities with labels a = 1 , 2 have 
asymptotic expansions around oo in the open half-planes C± (they are holomorphic expansions in the 
case of compactly supported measures da,dp). We will subsequently call the part of the expansion 
corresponding to negative powers of z or w, of a function f{z,w) the regular part of / and denote it 
{f(z,w))_ ^, {f(z,w))_ respectively. In all cases the regular parts are obtained by subtracting certain 
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polynomial expressions from functions holomorphic in C± and as such the regular parts arc holomorphic 
in these half-planes with vanishing limits at oo approach from within the respective half-planes. 

We will only indicate the main steps in computations for each entry, denoted below by (a, b). 
(0,1): With the help of the first approximation condition, we have 

q^(w)npo(z) 



ci[{w)Ilpo{z) = (^j . 



-my) 



w-y 

Using the Christoffel-Darboux Identities and the notation of Corollarv l4.3l we get 



qfHnpo(2)= (^J ■ 



{w + z){w — y) 



-my) 



{w + z){w — y) \J (w + z)(w ~ y) 

where we dropped the projection sign in the first term because A(—w) is a polynomial of degree one. 
Using now the partial fraction decomposition 

1 1/1 1 



{w + z){w — y) z + y \w — y w + zj^ 

we get that 

(q^H - ci^{ym-w)poi^)^^\ j {ali-z) - qg^(^))[n, {-z Y^)£]po(z) ^^^^^^ 



so 



{w + z){w-y) ) _^ \J {w + z){z + y) 

Observe that (-z- Y^)ipo(z) = 0, q^(-z)(-z- Y^)L and {y){-z-Y'^)L ^ -{y + z)q^{y)L 

(ciiw) - Ci^iy))Ai-w)p„iz)^\ ^( j (qg-(y))(z + Y^)Znpo(z) ^^^^\ 



{w + z){w-y) ' ) -w \j {w + z){z + y) 

/SiMllP£(f)d/^(,) = 0, 
J w + z 

because the /3 averages of q are zero. Thus 

{w + z)qf (w)npo(z) = c[^{w)A{~w)po{z). 
(2,0): Using the second Pade approximation condition and biorthogonality we easily obtain 

„T , ^TT , , R];^,{w)A{-w)po{z) + l 

R^„.(w)npo(z) = , 

w + z 

Now, substituting this formula into the formula for the third approximation condition, written as in 
equation (|5-5p . gives: 

„T , NTT ^^ R^.0(i«)A(-«;)po(z) - 1 

w + z 
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Restoring the collective notation of qa,Pa we obtain : 

{w + z)ci2{w)npo{z) = q^(w)A(-'U7)po(z) - 1. 

(0,1): To compute qj^(w)npi(z) we use the Padc approximation conditions (|5-6p . in particular the first 
condition gives us: 

q^{w)Ilpo{z)Wa{z) - q^{w)IlP^{z) - qj^HHR^ (z). 

We observe that this time we have to project on the negative powers of z. Thus the goal is to compute 
{q^{w)Upo{z)W^{z))_^. We have 

{w)Ilpo{z)da{x)\ _ ( f qo (w)A(-u;)po(z)da(x) ~ 



z — x ) -z \J {z — x)(w + z) 

q^{w)A{-w)po{x)da{x)\ ( /" q^(u;)A(-w)(po(z) - po(a;))da(a;) 



[z — x){w + z) /_r \J {z — x){w + z) 

We see that the first term is already regular in z. To treat the second term we perform the partial fraction 
expansion -, 7^ — \ — r = ^1 and observe that the term with does not contribute, while 

^ {z—x)i^w+z) w+x^z—x w+z^ z — x ' 

the second term 

q]^{w)A{-w){po{z) -po{x))da{x)\ _ f f q^ {w)A{-w){po{-w) - po{x))da{x) 



(w + x){'w + z) J-z yj {w + x){w + z) 

q^ {w)A{-w)po{x)da{x) 
{w + x){w + z) 

Thus 

q^{w)A{-w)L-'^Pi{z) q^{w)A{-w)L-'^pi{-w) 

Qo {w)i-ypi{z) = ■ ■ . 

w + z w + z 

In other words, 

{w + z)q^ {w)npi{z) = q^ {w)A{~w)L^^ {pi{z) - pi(-w)). 
More explicitly, the second term above can be rewritten as 

-qo (w)A(-w)Z"^Pi(-w) = qo (w)n j p{x)da[x). 

On the other hand 



cg{w)Ii j p{x)da{x)+ql{w)Ii jj 



Poix + y) ' J J Po{x + y) 

{w - y)p{x)da{x)df3{y) 



f3o{x + y) 



Now the second term qp (i(j)n J (■!" y)pix)da{x)d0(y) _ q ]3j,(,g^ygg q^(i(;)n(p(a;)|») is a projector on poly- 
nomials of degree < ?i — 1 and thus wqQ {w)Il{p{x)\l) — qQ {w)Il{p{x)\y) = w — w = 0, hence 

{w + z)q^ {w)Upi{z) = q^(w)A(-u;)pi(z), 
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where Pi{z) = L^^(pj(z) + ^(p|l)) as advertised earlier. 

(1, 1): We use again the Pade approximation conditions (|5-6p . this time muhiplying on the left by q[{w)Il 
and projecting on the negative powers of z , to obtain: 

With the help of the result for the (0, 1) entry, after carrying out the projection, we obtain 

{w + z)ciJiw)np^{z) = tiJ{w)Ai-w)p^iz)+ciJ{w)Ai-w) ( [ E^^MeI _ ^ (p|i; 

\J w + x Po 

We claim that 



w + x Pq 

Indeed, the left hand side of the equation equals: 

Po J J x + y Po J 1, 

^ I ^d/3(0 ^ -1. 

Po J w-^ 

Thus 

{w + z)q[ {w)Upi{z) = (w)A(-w)pi(2) - 1. 

(2, 1): This time we use projections in both variables, one at a time, and compare the results. First, let 
us use the projections in z. Thus 

(iz;)npi(z) = (q^Hnpo(2)W^a(z))__^ . 

Carrying out all the projections we obtain an expression of the form: 

r.T( mr. ( \ ^ qUw)A{~w)pi{z) Wc.{z) + F{w) 
H2 \W)^^Pi[z) — ; ; ■ 

w + z w + z 

Observe that, since q^iw) is 0{l/w) and the first term on the right is much smaller, F{w) = 0{1). More 
precisely, by comparing the terms at 1/w on both sides, we conclude that in fact, F[w) = 0{l/w). Now, 
we turn to the projection in w, resulting in an expression of the form: 

r.T( mr. f \ - qJWA(-w)pi(z) Wa*iw) + G{z) 

Vt2 \W)^^Pl\Z) — ; ; ■ 

w + z w + z 

This, and the fact that F{w) = 0{l/w), implies that F{w) = Wa-iw), G{z) = Wa{z). Hence 
(w + z)qf (u;)npi(z) = {w)A{~w)pi{z) - (W^„(z) + Wa'{w)). 
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(0,2): We use the projection in the z variable and the fact that by the Padc approximation condition 
(|5-5p . after exchanging a with /3, P2(2^) = ^\{z)Wf)-'{z) — Rq^. (z). Using the resuh for the (0, 1) entry 
we obtain: 

Ti \TT f _ ^oi'w)M~w)pi{z)Wf}.{z) f {w)A{-w)po{z)Wai3' {zY 

Qo [w)np2[z) — ■ ■ 

w + z \ w + z 

Carrying out the projection and reassembhng terms according to the definition of P2iz) we obtain: 

q^{w)A{-w)p2{z) q^{w)U{po\l) {w)A{-w)p2iz) 1 



Qo {w)Up2iz) 



w + z w + z w + z w + z 



(1,2): We use the projection in the z variable and the Padc approximation condition P2{z) = pi{z)Wp* (z)- 
Tiap'iz). Consequently, 

qf (uOnp2(z) = qJ{w)Upi{z)Wp,{z) - qj {w)URai3> (z) = 

qj {w)A{-w)pi{z) - 1\ /' T/ NTT / M^/ f W 
— — W^a.(,) - (qi {w)npo{z)Wai3'{z))_ ^. 

W ~r Z J ' 

Using the existing identities and carrying out the projection in the second term we obtain: 

[w + z)qf (i(;)np2(z) = qf (w)A(-u.)p2(z) - Wp^^z) - Wfiiw). 

(2, 2): The computation is similar to the one for (1, 2) entry; we use both projections. The projection in 
the z variable gives: 

T, ,TT , ^_ ^2k^)A{-w)p2{z) , F{w)-{W^.{w)+W^{z))Wfi.{z) + W^(i.{z) 
q^ \ w)vvp2yz) — 1 \ I . 

On the other hand, carrying out the projection in the w variable we obtain: 

q^(«;)np2(z) = q^H^(-^)P2(^) + G{z) - {Wp(w) + Wi,.{z))W^.{w) + Wp^.{w) ^ 
w + z w ~\- z 

Upon comparing the two expressions and using Lemma 15.21 we obtain F(w) ~ — H/q,-»^(u'), hence 

{w + z)q^(w;)np2(z) = q^(«;)A(-«;)p2(z) - W^.piw) - {W^, {w) + W^{z))W{^. (z) + M^„/3. {z) = 
q^(u;)[n, A(-u;)]p2 (z) - (M^„. (u;)W^. (z) + W^.fiiw) + Wf^.^^z)), 

where in the last step we used again Lemma 15.21 □ 
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